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Abstract. We derive explicit valuation formulae for an exotic path- 
dependent interest rate derivative, namely an option on the composition 
of LIBOR rates. The formulae are based on Fourier transform methods 
for option pricing. We consider two models for the evolution of interest 
rates: an HJM-type forward rate model and a LIBOR-type forward price 
model. Both models are driven by a time-inhomogeneous Levy process. 

1. Introduction 

The main aim of this paper is to derive simple and analytically tractable 
valuation formulae for an exotic path dependent interest rate derivative, 
namely an option on the composition of LIBOR rates. The formulae make 
use of Fourier transform techniques, see e.g. Eberlein, Glau, and Papapan- 
toleon (I2008p . and the change-of- numeraire technique. There are two models 
for the term structure of interest rates considered in this paper: a Heath- 
Jarrow-Morton-type forward rate model and a LIBOR-type forward price 
model, both driven by a general time-inhomogeneous Levy process. 

A standard approach to modeling the term structure of interest rates 



is that of Heath, Jarrow, and Morton (1992) In the Heath-J arrow-Morton 



(henceforth HJM) framework subject to modeling are instantaneous con- 
tinuously compounded forward rates which are driven by a d-dimensional 
Wiener process. However, data from bond markets do not support the use of 
the normal distribution. Empirical evidence for the non-Gaussianity of daily 



returns from bond market data can be found in Raible (2000, chapter 5) the 
fit of the normal inverse Gaussian distribution to the same data is particu- 
larly good, supporting the use of Levy processes for modeling interest rates. 
Similar evidence appears in the risk-neutral world, i.e. from caplet implied 



volatility smiles and surfaces; see Eberlein and Kluge (2006a 



The Levy forward rate model was developed in Eberlein and Raible (1999) 
and extended to time-inhomogeneous Levy processes in Eberlein, Jacod, 
and Raible (|2005jh In these models, forward rates are driven by a (time- 
inhomogeneous) Levy process; therefore, the model allows to accurately 
capture the empirical dynamics of interest rates, while it is still analytically 
tractable, so that closed form valuation formulae for liquid derivatives can 
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be derived. Valuation formulae for caps, floors, swaptions and range notes 
have been derived in Eberlein and Kluge (2006a, I2006bp . while estimation 



and calibration methods are discussed in 



Iberlein and Kluge| (|20U6at [20D7D . 



Moreover, Eberlein, Jacod, and Raible (2005) provide a complete classifi- 
cation of all equivalent martingale measures in the Levy forward rate model. 
They also prove that in certain situations - essentially, if the driving pro- 
cess is 1-dimensional - the set of equivalent martingale measures becomes a 
singleton. 

The main pitfall of the HJM framework is the assumption of continuously 
compounded rates, while in real markets interest accrues according to a dis- 
crete grid, the tenor structure. LIBOR market models, that is, arbitrage-free 
term structure models on a discrete tenor, were constructed in a series of arti- 



cles by Sandmann et al. (1995) , Miltersen et al. (1997) Brace et al. (1997) 



and Jamshidian (1997) In addition, LIBOR market models are consistent 
with the market practice of pricing caps and floors using Black's formula 
fcf. lBlack 19"76l) . 

Nevertheless, a familiar phenomenon appears: since the model is driven 
by a Brownian motion, it cannot be calibrat ed accurately to the whole te rm 
structure of volatility smiles. As a remedy, Eberlein and Ozkan (2005) de- 



veloped a LIBOR model driven by time inhomogeneous Levy processes. Val- 
uation met hods for caps and floors, usi ng approximation arguments, were 
presented in Eberlein and Ozkan (2005) and Kluge (2005) , while calibration 



issues for this model are discussed in Eberlein and Kluge (2007) 



The Levy forward price model is a market model based on the forward 
price - rather than the LIBOR rate - and driven by time inhomogeneous 
Levy processes; it was put forward by lEberlein and Ozkanl (2005, pp. 342- 
343). A detailed construction of the model is presented in |Kluge| (2005, 
Chapter 3); there, it is also shown how this model can be embedded in the 
Levy forward rate model. 

Although the forward LIBOR rate and the forward price differ only by an 
additive and a multiplicative constant, the two specifications lead to models 
with very different qualitative and quantitative behavior. In the LIBOR 
model, LIBOR rates change by an amount relative to their current level, 
while in the forward price model changes do not depend on the actual level 
(cf. |Kluge 2005 p. 60). There are authors who claim that models based on 
the forward process - also coined "arithmetic" or "Bachelier" LIBOR models 
- are able to better describe the dynamics of the market than (log-normal) 



LIBOR market models; see Henrard (2005) 



Another advantage of the forward price model is that the driving process 
remains a time-inhomogeneous Levy process under each forward measure, 
hence this model is particularly suitable for practical implementation. The 
downside is that negative LIBOR rates can occur, like in an HJM model. 

This paper is organized as follows: in section [2] we review some basic 
properties of the driving time-inhomogeneous Levy processes and in section 
[3] we describe the forward rate and forward price frameworks for modeling 
the term structure of interest rates. In section U] the payoff of the option on 
the composition is described and valuation formulae are derived in the two 
modeling frameworks. Finally, section [5] concludes. 
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2. TlME-INHOMOGENEOUS LEVY PROCESSES 

Let (O, J-, IF, IP) be a complete stochastic basis, where T = J~t* and 
the filtration F = (J r t)te[o.T*] satisfies the usual conditions; we assume that 
T* G M + is a finite time horizon. The driving process L = (£t)te[o,T*] is a 
time-inhomogeneous Levy process, or a process with independent increments 
and absolutely continuous characteristics, in the sequel abbreviated PIIAC. 
Therefore, L is an adapted, cadlag, real-valued stochastic process with in- 
dependent increments, starting from zero, where the law of Lt, t G [0,T*], 
is described by the characteristic function 
t 

B [e-] - e X p/ ( ib ,u - + /(e»- - ! - te ) M(ta ))d., (2.1) 

R 

where 6f G R, q G M+ and A^ is a Levy measure, i.e. it satisfies At({0}) = 
and J R (1 A |x| 2 )At(dx) < oo, for all t G [0, T*]. In addition, the process L 
satisfies Assumptions (AC) and (EM) given below. 

Assumption (AC). The triplets (t>t,ct,\t) satisfy 

T* 

J (\h\ +c t + J (1 A \x\ 2 )X t (dx)^j dt < oo. (2.2) 

R 

Assumption (EM). There exist constants M, e > such that for every 
u£ [-(l + e)M,(l + e)M] 

T* 

J J exp(ux)X t (dx)dt < oo. (2-3) 

{\ X \>1} 

Moreover, without loss of generality, we assume that JVu^n e^A^dx) < oo 
for all t G [0, T*] and all u G [-(1 + e)M, (1 + e)M\. 

These assumptions render the process L = (L t )o<t<T* a special semi- 
martingale, therefore it has the canonical decomposition (cf. Jacod and 
ShiryaevElSl II.2.38, and lEberlein et al. 2005|> 

t t t 

Lt = j b s ds + J ^/c~dW s + j j x(fi L -u)(ds,dx), (2.4) 

R 

where fi L is the random measure of jumps of the process L and W = 
(Wt)o<t<T* is a P-standard Brownian motion. The triplet of predictable or 
semimartingale characteristics of L with respect to the measure P, T(L\P) = 
(B,C,u), is 

t t t 

B t = f b s ds, C t = j c s ds, u([0,t] x^) = y j X s (dx)ds, (2.5) 

OA 

where A G B(R). The triplet (b, c, A) represents the local or differential 
characteristics of L. In addition, the triplet of semimartingale characteristics 
(B, C, v) determines the distribution of L. 
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We denote by 8 S the cumulant generating function (i.e. the logarithm 
of the moment generating function) associated with the infinitely divisible 
distribution with Levy triplet (b s ,c s , X s ), i.e. for z £ [— (1 + e)M, (1 + e)M] 

9 s {z) := b s z + ^z 2 + f (e zx - 1 - zx)X s (dx). (2.6) 

M 

Subject to Assumption (EM), 6 S is well defined and can be extended to the 
complex domain C, for z G C with $t.z € [— (1 + e)M, (1 + e)M] and the 
characteristic function of Lt can be written as 

f 

E [e iuLt ] = exp J 6 s (iu)ds. (2.7) 
o 

If L is a (time- homogeneous) Levy process, then (b s , c s ,X s ) and thus also 6 S 
do not depend on s, and 6 equals the cumulant generating function of L\. 

Lemma 2.1. Let L = (Lt)o<t<T* be a time-inhomogeneous Levy process 
satisfying assumption (EM) and f : M + — > C a continuous function such 
that |5ft(/)| < M. Then 



E 



Z -i I 

expJf(s)dL s =exp J 9 s (f(s))ds. (2.8) 



o 



(The integrals are to be understood componentwise for real and imaginary 
part.) 

Proof. The proof is similar to the proof of Lemma 3.1 in Eberlein and Raible 



[T9991 see also Kluge (2005, Proposition 1.9) □ 



3. Levy term structure models 

In this section we review two approaches to modeling the term structure 
of interest rates, where the driving process is a time-inhomogeneous Levy 
process. 

3.1. The Levy forward rate model. In the Levy forward rate framework 
for modeling the term structure of interest rates, the dynamics of forward 
rates are specified and the prices of zero coupon bonds are then deduced. 
Let T* be a fixed time horizon and assume that for every T G [0,T*], there 
exists a zero coupon bond maturing at T traded in the market; in addition, 
let U G [0,T*]. 

The forward rates are driven by a time-inhomogeneous Levy process 
L = (Lt)t£[Q t T*] on the stochastic basis F,P) with semimartingale 

characteristics {B,C,i>) or local characteristics (6, c, A). The dynamics of 
the instantaneous continuously compounded forward rates for T £ [0, T*] is 
given by 

t t 

f(t, T) = /(0, T) + J a(s, T)ds - J a(s, T)dL s , < t < T. (3.1) 



COMPOSITIONS IN LEVY TERM STRUCTURE MODELS 



5 



The initial values /(0, T) are deterministic, and bounded and measurable 
in T. In general, a and a are real-valued stochastic processes defined on 
SI x [0, T*\ x [0, T*\ that satisfy the following conditions: 

(Al): for s > T we have a(ui; s, T) = and cr(ui; s, T) = 0. 

(A2): (u, s, T) i y a(uj; s, T), a(uj] s, T) are V <g> B([0, T*])-measurable. 

(A3): S(u>) := sup SjT < T »(|a(u;; s, T)\ + \a(u; s,T)\) < oo. 

Then, ()3.ip is well defined and we can find a "joint" version of all f(t,T) 
such that (u;t,T) i-> /(t, T)(o;)l {t < T} is O <g> B([0, T*])-measurable. Here P 
and O denote the predictable and optional <7-fields on 17 x [0, T*\. 

Taking the dynamics of the forward rates as the starting point, explicit 
expressions for the dynamics of zero coupon bond prices and the money 



market account can be deduced (cf. Proposition 5.2 in Bjork et al. 1997). 



From Eberlein and Kluge (2006b, (2.6)), we get that the time-T price of a 



zero coupon bond maturing at time U is 



(T T \ 

J E(s, T, U)dL s - J A(s, T, U)ds J , (3.2) 
o o / 

where the following abbreviations are used: 

£(s,T,[/) :=S( S ,LT)-S( S ,T), 
A(s,T,U) :=A(s,U)-A(s,T), 

and 

T T 

A(s,T) := J a(s,u)du and S(s,T):= J a(s,u)du. 

sAT sAT 



Similarly, using Eberlein and Kluge (2006b, (2.5)) , we have for the money 
market account 

(T T \ 

J A(s,T)ds- J E(s,T)dL s . (3.3) 
/ 

In the sequel we will consider only deterministic volatility structures. 
Therefore, £ and A are assumed to be deterministic real-valued functions 
defined on A := {(s,T) £ [0, T*] x [0,T*];s < T}, whose paths are con- 
tinuously differentiable in the second variable. Moreover, they satisfy the 
following conditions. 

(Bl): The volatility structure £ is continuous in the first argument 
and bounded in the following way: for (s, T) £ A we have 

< £(s,T) < M, 

where M is the constant from Assumption (EM). Furthermore, we 
have that £(s, T) / for s < T and S(T, T) = for T G [0, T*]. 
(B2): The drift coefficients A(-,T) are given by 

A(s,T) = 6 s (Z(s,T)), (3.4) 

where 9 S is the cumulant generating function associated with the 
triplet (b s ,c s , A s ), s G [0, T}. 
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Remark 3.1. The drift condition (|3.4p guarantees that bond prices dis- 
counted by the money market account are martingales; hence, P is a mar- 
tingale measure. In addition, from Theorem 6.4 in Eberlein et al. (2005) , we 
know that the martingale measure is unique. 

3.2. The Levy forward price model. In the Levy forward price model 
the dynamics of forward prices, i.e. ratios of successive bond prices, are 
specified. Let = To < T\ < ■ ■ ■ < T/v < T/v+i = T* denote a discrete tenor 
structure where <5j = T + i — T, j G {0,1,... , TV}; the model is constructed via 
backward induction, hence we denote by T* := T/v+i-j for j G {0, 1, ... , N+ 
1} and 5* := S N+1 -j for j G {1, . . . , N + 1}. 

Consider a complete stochastic basis (Q, F, ) and let L = (Lt)t£[o,T*] 
be a time-inhomogeneous Levy process satisfying Assumption (EM). L has 
semimartingale characteristics (0, C, v T *) or local characteristics (0,c, X T *) 
and its canonical decomposition is 

t t 

L t = J y/c a &WT + f J X ^ L ~ v T l(ds,dx), (3.5) 

OR 

where W T * is a Py. -standard Brownian motion, fi L is the random measure 
associated with the jumps of L and v T * is the Py» -compensator of \i L . 
Moreover, we assume that the following conditions are in force. 

(FP1): For any maturity Tj there exists a bounded, continuous, deter- 
ministic function r](-, T) : [0, T] — ► M, which represents the volatility 
of the forward price process F(-, Tj, Ti+Si). Moreover, we require that 
the volatility structure satisfies 

i 

|£>(s,T fe )| <M, Vi€{l,...,N], 
k=l 

for all s G [0, T*], where M is the constant from Assumption (EM) 
and n(s, Ti) = for all s > Tj. 
(FP2): The initial term structure 5(0, T), 1 < i < TV + 1 is strictly 
positive. Consequently, the initial term structure of forward price 
processes is given, for 1 < i < N, by 

F(0,Tj,Tj + <5j)- m Tl) 



B(0,Ti + Si) 

The construction starts by postulating that the dynamics of the for- 
ward process with the longest maturity F(-, Tj* , T*) are driven by the time- 
inhomogeneous Levy process L, and evolve as a martingale under the termi- 
nal forward measure Py* . Then, the dynamics of the forward processes for 
the preceding maturities are constructed by backward induction; therefore, 
they are driven by the same process L and evolve as martingales under their 
associated forward measures. 

Let us denote by Pt* the forward measure associated with the settle- 
ment date T*_ 1 , j G {!,..., N + 1}. The dynamics of the forward price 
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process F(-,T?,T*_ 1 ) is given by 

\o 
where 

t t 
L^- 1 = J ^dWp- 1 + J J x(/i L -v T j-i)(ds,dx) 

OR 

is a time-inhomogeneous Levy process. Here W is a P7 1 * 1 -standard 
Brownian motion and v 3-1 is the Py* i -compensator of ji . The forward 
price process evolves as a martingale under its corresponding forward mea- 
sure, hence, we specify the drift of the forward price process to be 

b(s,T*,T*_ l ) = -l(r l (s,T*)fc s 

' e ri(s,T*)x _ 1 _ v ( 8t Tf)x) X?-\dx). (3.6) 



The forward measure Pr*_ 1 , which is defined on (0, J 7 , (^i)o<t<Tf_ 1 )i is 
related to the terminal forward measure Py* via 

dP T * F^T*,^) BfaTf^U ^i-i' fc ' 

In addition, the Pt» -Brownian motion is related to the Pr* -Brownian 
motion via 

i 





t 



W - >>( s ,r fe *) K^d s . (3.7) 





Similarly, the Py* i -compensator of fi L is related to the Py* -compensator 
of /i L via 

v^i~ x (ds, dx) = exp ^77(s,T J *_ 1 )x^i/ T J- 2 (ds,da;) = ... 

= exp (x^r)(s,TZ)\ v T *(ds,dx). (3.8) 

Remark 3.2. The process L = L T * , driving the most distant forward 
price, and L^'- 1 , driving the forward price F(-,Tj,Tj_^), are both time- 
inhomogeneous Levy processes, sharing the same martingale parts and dif- 
fering only in the finite variation parts. Applying Girsanov's theorem for 
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semimartingales yields that the Py* i -finite variation part of L is 
J c s 2^ ds + j Jx(e W (x^77(s,2fc))-l] ^ T *(d S ,dx). 



\k-l / OR 

4. Valuation of options on compositions 

Consider a discrete tenor structure = To < T\ < ■ ■ ■ < T/v < Tjy+i = 
T*, where the accrual factor for the time period [T,T + i] is 5i = 2i+i — T, 
i 6 {0, 1, . . . , N} and let L(sj, T) denote the time-Sj forward LIBOR for the 
time period [Tj,T + i]. The composition pays a floating rate, typically the 
LIBOR, compounded on several consecutive dates. The rates are fixed at 
the dates Si < T and the value of the composition is 

N 

H(l + 5 i L(s i ,T l ))- 

i=l 

therefore, the composition equals an investment of one currency unit at the 
LIBOR rate for N consecutive periods. The value of the composition is 
subjected to a cap (or floor) denoted by K and is settled in arrears, at time 
T*. Hence, a cap on the composition pays off at maturity the excess of the 
composition over K, i.e. 

+ 

'l + 6iL(8i,T$-K\ , 
and similarly, the payoff of a floor on the composition is 

Notice that without the cap (resp. floor), the payoff of the composition would 
simply be that of a floating rate note, where the proceeds are reinvested. 
Similarly, if we only consider a single compounding date, then we are dealing 
with a cap let (resp. floorlet), with strike % := — f^. 

In the following sections, we present methods for the valuation of a cap on 
the composition in the Levy-driven forward rate and forward price frame- 
works. The value of a floor on the composition can either be deduced via 
analogous valuation formulae or via the cap-floor parity for compositions, 
which reads 

C(T*;K) =F(T*;K) + B(0,T 1 ) - KB(0,T*). 

Here C(T*;K) and F(T*;K) denote the time- To value of a cap, resp. floor, 
on the composition with cap, resp. floor, equal to K. 

4.1. Forward rate framework. In this section we derive an explicit for- 
mula for the valuation of a cap on the composition in the Levy forward rate 




model, making use of the methods developed in Eberlein et al. (2008) As a 
special case, we get valuation formulae for caplets in the Levy forward rate 
framework that generalize the results of Eberlein and Kluge (2006a) since 
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we do not require the existence of a Lebesgue density (which is essential in 
the convolution representation of option prices; cf. Raiblc 2000, Chapter 3). 

Firstly, we calculate the quantity that appears in the composition. By an 
elementary calculation, we have that 



B( Si ,Ti) B(0,Ti\ 



cxp 



A(s,Ti,T i+1 )ds - f E(s,Ti,T i+1 )dL s 
o 



B{s u T l+1 ) B{0,T i+l 
Using the fact that 1 + 5iL(si,Ti) = R ?f s '^ T '' > ; we immediately get 



i=l 



5(0, T*) 

N 



(n y n y 

/ A(s,Ti,T i+1 )ds-J2 / S(s,T i5 r i+1 )dL 
1=1 n i=l n 



Next, we define the forward measure associated with the date T* via the 
Radon-Nikodym derivative 



dP r . 
dP 



1 



B$lB(0,T*) 



(T* T* \ 

-J A(s,T*)ds + J S(s,T*)dL s J . 



The measures P and Py* are equivalent, since the density is strictly positive; 
moreover, we immediately note that E[ gM ^ Q ^y ] = 1. The density process 

related to the change of measure is given by the restriction of the Radon- 
Nikodym derivative to the cr-field Ft, t < T* , therefore 



E 



dP r * 



dP 



B(t, T* 



BfB(0,T*) 



-J A(s,T*)ds + J S(s,T*)dL s . 



This allows us to determine the tuple of functions that characterize the 
process L under this change of measure and we can conclude, using The- 
orems III. 3. 24 and II. 4. 15 in Jacod and Shiryaev (2003) , that the driving 
process L = (£t)te[o,r*] remains a time-inhomogeneous Levy process under 
the measure Pj 1 * . 

According to the first fundamental theorem of asset pricing the price of an 
option on the composition is equal to its discounted expected payoff under 
the martingale measure. Combined with the forward measure defined above, 
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this gives 



C(T*;K) = E P 



will 



B( Si ,Ti) 
B™ \^B{s h T i+1 ) 

N 



K 



£(0,T*)E Pt » 



n 



B( Si ,Ti) 



K 



±B(s u T i+1 ) 

= £?(0,T*)E Pt , [(e W H-K) + ] , 
where the random variable H is defined as 

N s ' N s * 

H:=log +^ J A(s,T u T l+1 )ds-Y^ / T h T i+1 )dL s . 

i=l i=l 

Let us denote by the moment generating function of H under the 
measure Py. The next theorem provides an analytical expression for the 
value of a cap on the composition. Preceding that, we provide an expression 
for M^j (z) for suitable complex arguments z. 

Lemma 4.1. Let M and e be suitably chosen such that S(s, T) < M' for 
all s,T G [0,T*\ and E{s,T i+1 )l [s . jS . +l] (s) < M" for all s, Si ,T i+1 G [0,T*], 
where < M" < M' < M and Mr > N + 1. T/ien, for each R G J 2 = 



[1 



A/-A/"(Af+l) i _,_ M-M"(N+1) 
M'+M"{N+l) ■ 



,1 + T7T- 



or 



every z G C wi/i = i? 



M T U * (z) = Z z exp U (^z Ms, T it T i+1 )l [0>Si] (s) - 8 S (e( S , T 

+6 S (Z(s, T*)-zJ2 r i} T i+1 )l [0jS .]( a )) j ds J 

wnere . — ■ 

Proof. Fix an i? G I2. Then, for z G C with = R, and denoting by 

M s , T ) = EiIo s ( s ' T i+i) 1 h,^+i]( s )' we § et that 

N 

m(-zJ2 t u r i+1 )i [0iSi] ( s )) + s( s , r*) 

i=l 

N 

= M(z^X(s,T i+1 )l [suSi+l] (s)-zX(s,T*)) +S( S ,T*) 



i=0 



= »^(l-z)(E( Sj r*)-E( S ,T))J + £>,T) 
< \l-R\\Z(s,T*)-Us,T)\ + \Us,T)\ 

* m'Im^nVi) {M ' + M " {N + 1)} + A///(iV + 1} = AL (41) 
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Now, define the constants 



Z := exp 2 lo, 



fl(0,?i) 
B(0,T*) 



f N 

n i=l 



and i?i := iJo x exp ^— J o r ^4(s, T*)ds^ . Hence, the moment generating 
function of H is 

Mj£*(z) \exp(zH) 



Ep 7 



exp (^( log §jb^ + £ / A ( a > T *> T *») d 

N s » 

Z(s,Ti,T t+1 )dL 



exp L- J A(s,T*)d}j x 

/ N Si 

exp -z^ J S( Sj T i ,T m )dL s + y £(s,T*)dL s 



x E P 



^-^^ s ( s > r i,rj + i)l[ 0)S .](s) + £(s,T*)^ dL s 

= Ziexp / I S ( - z^2H(s,T h T i+1 )l[ 0iSi ](s) + £(s,T*)J J ds, 
o ^ i=l / 

where for the last equality we have applied Lemma [2. 11 which is justified by 
flUD - In addition, we get that M^* (R) < oo for i? G T 2 . □ 

Theorem 4.2. Assume that forward rates are modeled according to the Levy 
forward rate model. The price of a cap on the composition is 



C(T*;K) 



2tt 



Mjj (R - iu) 



j^l+iu—R 



(iu - R)(l + iu-R) 



du, 



where Mjj* is given by Lemma \J7J\ and R £ (1, 1 + m'+M"(N+1) 1 ' 

Proof. Firstly, let us recall that the Fourier transform of the payoff function 
f(x) = (e x — K) + , K G R_|_, corresponding to a call option is 



/(*) 



K 



l+iz 



iz(l + iz) 



(4.2) 



for z£C with Qz G (1, oo) =: X±; cf. Example 3.15 in Papapantoleon (2007) 
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Now, since the prerequisites of Theorem 2.2 in Eberlein et al. (2008) are 



satisfied for R G 1\ fl Ti , we immediately have that 



C(T*;K) = B(0,r)Ep T , 
B(0, T* 



2tt 



Mfj (R - iu) 



j^l+iu—R 



{iu - R)(l + iu-R) 



du, 



and the assertion is proved. 



□ 



4.2. Forward price framework. The aim of this section is to derive an 
explicit formula for the valuation of a cap on the composition in the Levy 
forward price model. Once again, the valuation formulae will be based on 



the methods developed in Eberlein et al. (2008) 



We begin by noticing that the quantity that appears in the composition 
can be expressed in terms of forward prices, since 

l + 5 i L(;T i )=F(;T i ,T i+1 ), 

and the forward prices are the modeling object in this framework. We know 
that each forward price process evolves as a martingale under its correspond- 
ing forward measure; moreover, we know that all forward price processes are 
driven by the same time-inhomogeneous Levy process (see also Remark [32]). 
Therefore, we will carry out the following program to arrive at the valuation 
formulae: 

(1) lift all forward price processes from their forward measure to the 
terminal forward measure; 

(2) calculate the product of the composition factors; 

(3) price the cap on the composition as a call option on this product. 

Appealing to the structure of the forward price process and the connection 
between the Brownian motions and the compensators under the different 
measures, cf. equations (|3.7|) and (|3.8|) . we get that 

(t t 

/ t t 

= f(o, t;, r;_i) ex P J b{ s , r;, r*)ds + J v ( s , r/)dLf 



(4.3) 



for all j € {1, . . . , N + 1}. Here L T * is the driving time-inhomogeneous Levy 
process with Py* -canonical decomposition 

t t 
LT = JyfcdWT + J J x{» L ~ v T *)(ds,dx), (4.4) 
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and the drift term of the forward process Tj", Tj ti _ 1 ) under the terminal 
measure Pj-« , is 

b(s, T*, T*) = -c. fen{8, T*)) 2 + V (s, 27) £ ^, 

- 1 (( e ^ T /) - lyx&vfrTt) _ xr? ( s , T ;)) Af(dx). 

R 

(4.5) 

It is immediately obvious from (g3J), flO)) and g3]) that F(-,T*,T*^) is 
not a Py« -martingale, unless j = 1 (where we use the convention that 

EU = °)- 

Now, the composition takes the following form 

N N 

J] (1 + <5iL(si,Tj)) = J] F(s*,T*,T*_ 1 ) 
i=i 3=1 



B(0,T* N ) 



(4.6) 



x exp 



5(0, T* 

' N 3 N 3 1 

2 / 6( S ,T/,T*)d S + ^ / r?( S ,T;)dLf 
V' =1 J' =1 / 

where s* = sjy+i-j, j 6 {1, • • • , N}. Define the random variable 

N S Z N Sj 



H:=l °Z^7P^ + lt J Hs,T*,T*)ds + f2 J V(s,T*)dLT (4.7) 
j=i i =l 



and now we can express the option on the composition as an option depend- 
ing on this random variable. The next theorem provides a formula for the 
valuation of a cap on the composition. 

Theorem 4.3. Let forward prices be modeled according to the Levy forward 
process framework. Then, the price of a cap on the composition is 

B(0 T*) f xi+iu-R 

C(T*;K) = '- / M H (R - iu)- — — du, (4.8) 

v ' ' 2vr J K (iu — R)(l + iu — R) v ; 



where the moment generating function of H is given by Lemma \4-4\ an d 
*€(1,#]. 

Proof. The option on the composition is priced under the terminal forward 
martingale measure Pp*. Using (|4.6[) and (|4.7p . we can express the cap on 
the composition as a call option depending on the random variable H. Then 
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we get 



C{T*;K) = £(0,T*)Ep 3 



A 1 



U' =1 



S(0,T*)E Pt , [(e H -K)~ 
5(0, T* 



2tt 



M ff (i? - iu)- 



j^l+iu—R 



iu - R)(l +iu- R) 



du 



where we have applied Theorem 2.2 in Eberlein et al. (2008) and used (|4.2p 
once again. □ 

Lemma 4.4. Let M and e be suitably chosen such that | X^fcLi v( s i < 
M' for some M' < M and for all s G [0, T*} . Then, for each R G [0, ^] we 
/io«e i/iai Mjj(R) < oo, and /or ewer?/ z G C u>ii/i Siz G [0, jp] the moment 
generating function of H is 



M H (z) = Z z exp 



J ( z E b ^ T h T *)+ C (* £ vis, 27)) ) d, j , 



where Z = ^(q'j^j and #J* *s cumulant generating function associated 
with the triplet (0, c s ,Aj ). 



Proof. Fix an i? G [0, -^p-] and then for z G C with = R we get 



A? 



fe=i 



i? 



A 



k=l 



~ M' 



(4.9) 



Now, define the constant 



Zo 



B(0,T* 



I 



cxp 



Then, similarly to the proof of Lemma 14.11 we get 



N jt 

Y, J b(s,T*,T*)ds 



M H {z) = Z 2 TE Pt , 



T* 



exp / z^2<r)(s,T-)l [0iS *](s)dLT* 
Vo i =1 

= Z 2 exp (f 9r (zJ2 V (s,T*)y)j , 

where for the last equality we have applied Lemma [2. 11 which is justified by 
(|4.9p . Note also that rj(s, T*) = for s > s*, which is the fixing date for the 
rate; accordingly, b(s,T*,T*) = for s > s*, cf. (|4"3j) . 

In addition, we get that M H (R) < oo for R G [0, Mr]. □ 
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Figure 1 . Implied normal volatilities for a 5Y option on the 
3M LIBOR. 



4.3. Numerical illustration. In order to get an idea about the difference 
in the prices of an option on the composition of LIBOR rates between the 
classical, Brownian-driven, HJM model and the Levy forward rate model, 
we set up an artificial but reasonable market. Our aim here is not to give 
a complete analysis but rather a flavor of the impact we can expect. As an 
example, we look at a 5 year caplet on a composition of LIBOR rates. 

Rates are assumed to be flat at 4% across all maturities and instruments. 
Moreover, we assume that all volatilities are marked according to the SABR 
model (cf. |Hagan et al. 2002[ ) with parameters a = 1%, a = 40%, (3 = 0% 
and p = 30%. For the calibration, we choose the Vasicek volatility struc- 
ture with a fixed parameter a = 0.05. Not surprisingly, the normal inverse 
Gaussian Levy model calibrates well to the market smile whereas the clas- 
sical HJM model only gets the ATM point right but cannot reproduce any 
smile or skew; see Figure [TJ 

The results for a 5Y caplet on a composition of LIBOR rates are shown 
in Figure El It should not be surprising that the classical model does not 
produce any smile whereas the Levy model does; note that the ATM prices 
are also different, where ATM ~ 1.22. Moreover, as has been observed in 
many other situations, the Brownian-driven model overprices the ATM op- 
tions and underprices the in- and out-of-the-money options compared to the 
more realistic Levy model. 



5. Conclusion 

We have presented valuation formulas, based on Fourier transforms, for 
pricing an option on the composition of LIBOR rates in the forward rate and 
forward price models driven by time-inhomogeneous Levy processes. Anal- 
ogous formulas can also be derived for the affine 'forward price'-type frame- 
work proposed by Keller-Ressel, Papapantoleon, and Teichmann (2009) ; this 
framework combines the analytical tractability of the forward price frame- 
work described here with positive LIBOR rates. 
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1.45 



Figure 2. Implied normal volatilities for a 5Y option on composition. 



The challenge ahead is to derive valuation formulas in the LIBOR model 
driven by a time-inhomogeneous Levy process. This task requires some so- 
phisticated approximations due to the structure of the dynamics in LIBOR 
market mo dels; the interested reader is referred to|Siopacha and Teichmann 
(|2007p and Hubalek, Papapantoleon, and Siopacha (2009) for a detailed anal 
ysis. 
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